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Abstract. In the paper rue construct an hierarchy of integrable 
Hamiltonian systems inhich describe the variation of n-iuave en- 
velopes in nonlinear dielectric medium. The exact solutions for 
some special Hamiltonians are given in terms of elliptic functions 
of the first kind. 



1. Introduction 

The describtion of the nonlinear n-uiave interaction plays important 
role in many areas of physics including optics, accoustic, plasma and 
fluid physics, see e.g. |Hol08l IDH92I IALMR98I IKRB79] . For example 
one can find the solutions of three mave equations in early ruorks in 
the area of nonlinear optics [ABDP62J. 

In this paper rue study a hierarchy of integrable Hamiltonian sys 
terns on the space of linear maps L('H_,'H + ) betrueen tino complex 
finite dimensional Hilbert spaces "H_ and T-L + ruhich is the particular 
case of a more general hierarchy defined on the Banach Lie Poisson 
space related to the restricted Grassmannian, see |GOlO| . This hier 
archy, as rue mill shorn, can be used to describtion of the variation of 
plenary ujave envelopes through the dielectric nonlinear medium. 

In particular rue study in detail the case ruhen dimH + = 2 and 
dim"H_ = 3. The case ruhen dimH + = 2 and dimH- = 2 is also 
investigated up to giving explicite formulas for solutions. Hamiltoni 
ans in these truo cases describe a nonlinear optical system consisting 
of six and four ruaves interacting in a nonlinear medium displaying 
Kerr like effects and causing convertion of light modes. 

For the paper self sufficiency rue give in Section [2] a short presen- 
tation of the ruave optics background concerning the sloru moving 
ruave envelopes in the nonlinear medium. 

l 
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In Section [3] the construction of the hierarchy of integrable Hamil- 
tonian systems mentioned above is given and is shoiun that it has a 
rich family of integrals of motion in involution. 

In Section[4]uje investigate the (2+3)-dimensionaI system. Applying 
reduction procedure me construct the angle-action coordinates for this 
systems and integrate it in quadratures. 

The (2 + 2) dimensional case rue consider in Section [5] We shorn 
that this case can be obtained as a reduction of the previous one and 
find solutions for it in terms of elliptic functions of the first kind. 

To conclude the paper me shortly discuss in Section [6] the physical 
interpretation of separate terms of one of the Hamiltonian in (2 + 3)- 
dimentional case. 

2. Wave propagation in dielectric media 

The nonlinear mave optics deals mith the interaction of the eletro- 
magnetic maves mith the medium. The complicated character of this 
interaction is manifested in the nonlinear dependence of the medium 
polarization field P(t, x) on the electric field E(i, x). Since the efects 
caused by the magnetic field are much meaker then the ones for mhich 
the electric field is responsible, so, they are usually neglected in mave 
optics problems, see e.g. [BC90J. 

Let us start from Maxmell equations in dielectric medium mithout 
the conduction current 

VxE(t) = ~B(t) (2.1) 
V-(e„E(t) + P(t)) = (2.2) 

V x H(t) = e ^E(t) + J^P(t) (2.3) 

V • B(t) = 0, (2.4) 

inhere ^-(t) is the polarization current. In optics one usually has deal 
mith nonmagnetic media mhat is expressed in the dependence 

B(t) = // H(t) (2.5) 

betmeen the magnetic induction and the magnetic field. On the other 
hand the polarization P(t) dependence on the electric field has in 
general the nonlinear functional character 

P(i)=P[E(i)]. (2.6) 

Substituting (|2.5|l into Q2. 1 ^ 2.4 ) one obtains 

V x (V x E(t)) = ~~m - ^o|p(t). (2.7) 
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Next, expressing E(£) and P(t) in terms of their Fourier transforms 
E(o>) and P(w) one reiurites equation (|2.7p and (|2.6p as folloms 



V x V x EM 



E(w) + u;VoP(u/) 



W ~* — (2.8) 

(2.9) 

,E(w„,)), 
(2.10) 

e.g. see [BC90J, mhere the n linear form x (^l, ■ ■ ■ i^n), called n th - 
susceptibility tensor, describes the order of nonlinearity of the inter- 
action of the electric field mitti the medium. 

The ones of possible solutions of (|2.8p - (|2.9p are the running plane 

maves 

E(w) = E(w, n)e iK (2.11) 



PH = P[%)]. 

The equation (|2.9p in the most general setting has form 

oo 

P(w) = ^2e dtd! ... dcu n x {n) (^i,---,^n)(E(a;i), 



inhere k := k ■ r and k ■ E(o>, k 
obtain 

d 2 - , , _ . d 



0. Substituting (T2~TT) into ((278) uje 

1 



-E(w, «) + 2i—E{u, k) = E + T^P[E(a;, «;)e* K ]e 



IKl^—lK 



(2.12) 



If the variations of luave envelopes E(w, «) are sufficiently sIouj in 
variable k then me can assume 



(2.13) 



and neglect the second derivative in (|2.12p . 

So, in the slomly varying approximation me have a system of first 
order differential equations 











« 






an 



(2.14) 



on the infinite family of functions E(w, k) parametrized by the real 
numbers wGl 

In practical applications one restricts usually its attention to the fi- 
nite number of the running plane maves mhich are labeled by the fre 
quencies u)i, . . . , and polarizations of the envelopes E(wi, k),. . .,E(un, k). 
In this case the systsem of equations (|2.14p reduces to the finite system 
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of Hamilton equations 



d .dH 



(2.15) 



dn dzi 
d _ .dH 



(2.16) 
(2.17) 



dn 1 dzi 



mith the Hamiltonian H = H(zi, . . . , Z2N, ■ ■ ■ %2n) mhich depends 
in nonlinear may on the independent modes Z\, . . . , Z2N describing a 
system of N running plane maves mhich slomly vary in space and the 
evolution of a system of nonlinearly interacting harmonic oscillators 



In the next section me mill construct a hierarchy of such type inte 
grable Hamiltonian systems. 

3. Hierarchy of Hamiltonian integrable systems on Matj\/ X Ar(C) 

In the paper |GOl0j me have constructed a hierarchy of Hamilton 
ian integrable systems on L 2 (H-,'H+), mhere "H_ and T-L + are com- 
plex separable Hilbert spaces (finite or infinite dimensional) and L 2 
denotes the class of Hilbert Schmidt operators. In this paper me mill 
restrict our considerations to the subcase mhen both fi_ and H + are 
finite dimensional. In that case many analytical difficulties vanish and 
L 2 (H,_,H, + ) coincides mith the set of all linear operators L(7i-,1-L+). 
In subsequel me mill identify L{U-,H+) = Mat MxN (C) = C MN mith 
the M x N complex matrices or mith <C MN mhere dime H- = M and 
dime H+ = N. 

We define the Poisson bracket on the space of smooth functions 

C°°(M at MxN (*£■)) in the standard may, i.e. 



mhere Z G MatMxJv(C) and its conjugate Z + G MatNxM(C). 

In order to obtain a family of Hamiltonians in involution me need 
to introduce the folloming notation. Let TL := H+ ©"H and P± be the 
orthogonal projector on 1-L±. We mill consider a self adjoint operator 
// G L(T-[) in the block form consistent mith the above decomposition 
of Ti as 



mhere operators A G L(7i + ) and D G L("H_) are self-adjoint. More- 
over me mill assume that eigenvalues of both A and D differ from 



Kum90j. 




(3.1) 




(3.2) 
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each other. In the follouuing A and D mill be interpreted as parame 
ters of considered Hamiltonian hierarchy. 

Proposition 3.1. The functions 

H ktX := Tr(/i + AP+) fc , (3.3) 
where A G 1 and k E N, are in involution 

{Hk,\, Hn,\>} = (3.4) 
with respect to the Poisson bracket (|3. \) . 

Proof. Without loss of generality me can assume that A' = . Since 
the derivative DHk,\{p) of H^x is given by 

DH kt M = k(n + \P + ) k - 1 (3.5) 

me gather that 

ftT-f 

- ikP_{fi + \P + ) k - 1 P + . (3.6) 



dZ 

Thus me obtain 

{H k ^H nfi } = knTi (P-(ji+ \P + ) k - l P + ^ n - l P-- (3.7) 

-P_/i"- 1 P + (/i + \P + f- l P_) = 
= iknTr ((// + AP+)*" 1 P+f/ 1 - 1 - /i n-1 P+(/i + XP+) k ~ 1 )- 

-iknTr (P + (fi + \P + ) k - l P + iJi n - 1 - P+^P+Qi + AP+) fc - 1 ). 

Note that due to the invariance of Tr under cyclic permutations the 
second term in the expression above vanishes identically. Thus me 
get 

{H ktX ,H nfi } = iknTr( f i n - 1 [( f i + \P + ) k - 1 ,P + }) = (3.8) 

= UhnTr (M n_1 [(At + AP+)*" 1 , fi + AP+]) - 
A 

-jiknTr (/i^K/i + AP+) fe -\ //]) = 
A 

= --i&nTr ((/i+AP + ) fc - 1 [/ i ,^- 1 ]) = 0. 
A 

□ 

Let us note here that functions -P^.a defined in (|3.3|l form an in 
tegrable hamiltonian hierarchy. This hierarchy possesses tmo extra 
families of integrals of motion, mhich are analogues of Manley Rome 
integrals, see e.g. [H0IO8 . 
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Proposition 3.2. The functions 

a k := Tr(A k ZZ + ) 5 k := Tr(D k Z + Z), (3.9) 
k e N, commute with each other and with Hamiltonians (|3.3|| 

{a k , ai} = {a k , 5 t } = {5 k , S t } = 0. (3.10) 
Proof. We note that 

^ = Z+A k ^ = A k Z (3.11) 
dZ 8Z+ v ; 

^k =D k z+ ^L = ZD k (312) 
dZ dZ+ v ' 

Thus iue have 

{a k , aft = i Tr(Z + A k A l Z - A k ZZ + A l ) = (3.13) 
{5 k ,5i} = iTi(D k Z + ZD l - ZD k D l Z + ) = (3.14) 
{a k ,6i} = iTi(Z + A k ZD l - A k ZD l Z + ) = 0. (3.15) 
Moreover rue have 

{a k ,Hi >x } = UTr (Z+ A k P+(/i + AP+)'- 1 P_A k 'ZP-fji + XP^ 1 P+) = 

(3.16) 

= il Tr {P^(P + ^P + ) k (n + XP+Y^ 1 - (P + /iP+) V-O* + AP+y- 1 ) = 
= 2/Tr (fi(P +f iP + ) k (fi + AP+)'- 1 - (P + /iP + )V(/^ + AP+)^ 1 )- 

-z/Tr (p + fi(p +f iP + ) k (fi + xp+y- 1 - (p + fiP + ) k fiP + (fi + xp+y- 1 ) = 

= il Tr (/,[(P +/ iP + ) fe , (/i + AP+y- 1 ]) - il Tr ((P +/ ,P + ) fc+1 (/i + AP+)' -1 - 

-{P + ^P + ) k+ \^XP + ) 1 - 1 ) = */Tr {(n+XP + )[(P +f iP + ) k , ( /i+ AP + )'- 1 ])- 

-UX Tr (P + [(P +/ ,P + ) fe , (/x + AP+y- 1 ]) = 0. 

Proof that {5 k , H^\} = is analogous to (|3.16p . □ 

The Hamilton equations generated by any H G C 00 (MatMxN{ < &)) 
given by the Poisson bracket (|3.1|l are the follovuing 

• OH • , OH . 

z = 'gz^ z =~<8z- (317 » 

In the case ujhen H = H ki \ they assume the form 

Z = kP + {n + AP + ) fc " 1 P_ (3.18) 

Z + = kP_{fi + XP + ) k - 1 P + . 

These equations are in general nonlinear. Even if the family of in 
tegrals of motion in involution is big enough, it may be technically 
difficult to find their explicit solutions. Thus in next sections we mill 
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restrict our considerations to several specific cases uuhen it is possible 
to solve the system in the explicit may. 

Combining the Hamiltonians (|3.3p for A = and k = 4, 5 

# 4)0 = Tr /i 4 = Tr (A 4 + D A + AD 2 Z + Z + AAZDZ + + (3.19) 

+ AA 2 Z Z + + 2(Z + Zf) 

H 5fi = Tr ^ = Tr (A 5 + D 5 + hD 2 Z + AZ + hDZ + A 2 Z+ (3.20) 

+ 5A(ZZ + ) 2 + hD(Z + Z) 2 + 5A 3 ZZ + + 5D 3 Z + Z) 

mith Hamiltonians otk and 5k for A; = 1,2,3 me find that the tmo 
folloming functions 

H = - Tr(Z + Z) 2 + Ty{AZDZ + ) (3.21) 
F = Tr A(ZZ+) 2 + Tr(D 2 Z+AZ + DZ + A 2 Z) (3.22) 

belong to the hierarchy of Hamiltonians generated by (|3.3|l and (|3.9p . 

Equation of motion mith respect to Hamiltonian \?>.2\) are the fol 
loming 

Z = AZD + ZZ + Z (3.23) 

Z + = DZ + y4 + Z + ZZ + . 

They can be considered as a pair of coupled Ricatti type equations 
on the variables Z and Z + . 



4. Solution in (2 + 3) dimensional case 

In this section me consider in details the Hamilton equations in the 
case mhen dim'H + = 2 and dim"H_ = 3. Without loss of generality 
me can assume that the matrices A and D are diagonal. So the matrix 
H defined in (|3.2p assumes the form 



( o-i 





Zl 


z 2 







a 2 


Vl 




^3 


Z\ 


Vi 


d\ 








Z2 


V2 





(h 





\Z3 


V 3 








d 3 ) 



(4.1) 



mhere Zk, Vf. G C and a^,dk G R. 
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Expressing the integrals of motion H, F, a 0l a \, 8\, 8 2 m the vector 
coordinates z = (zi,z 2 ,z 3 ) T and v = (vi,v 2 ,v 3 ) T we obtain the fol- 
lowing six functionally independent integrals of motion in involution: 

H = \{ z+z ) 2 + \{ v+v f + \ v+z \ 2 + a i z+Dz + d2V + Dv (4.2) 
F = ai(z + z) 2 + a 2 (v + v) 2 + (di + a 2 ) \v + z\ 2 + z + zz + Dz + v + vv + Dv+ 

+ v + zz + Dv + z + vv + Dz + a\z + Dz + a\v + Dv + a 1 z + D 2 z + a 2 v + D 2 v 

(4.3) 

«o = <^o = z+ z + v + v (4.4) 
a\ = a\z + z + a 2 v + v (4.5) 
81 = z + Dz + v + Dv (4.6) 
S 2 = Z +D 2 z + v + D 2 v. (4.7) 

Nom me define for k = 1, 2, 3 functions 

r] k :=v k z k , r k := \z k \ 2 - \v k \ 2 , s k := \z k \ 2 + \v k \ 2 (4.8) 

on the phase space C 6 . They span the Lie algebra mith respect to the 
Poisson bracket l|3.1jl : 

{Vk,m} = ir k 5 kl 
{r k ,r]i} = 2irj k 8 k i 

{r k ,fji} = -2ifj k 5 H (4.9) 

{s k , si} = {s k , n} = {s k , rji} = {s k , fji} = 

{Vk,Vi} = {Vk,Vi} = 

mhere k, I = 1, 2, 3. 

Let us consider Lie Poisson space u(2)* dual to the Lie algebra 

u(2) := {X G Mat 2x2 (C) \ X + + X = 0} 

of the unitary group U(2). We identify u(2)* mith u(2) by the pairing 

(X ; Y) := Tr(XY) (4.10) 

and introduce the folloming coordinates 

X = i(j (4.11) 

on u(2)*. In these coordinates the Lie Poisson bracket 

{f,g} LP (X) := (X ; [Df(X),Dg(X)\) (4.12) 
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assumes the follouuing form 



{/■«M«.™*=Mr ('f!-H') + (4.13) 



dr drj drj dr J \ drj dr dr drj 
Note that functions Si, s 2 , s 3 , and 

c(r,r),fj) := — + 2 \r]\ 2 (4.14) 

are Casimirs for the Lie Poisson bracket (|4.13|l . 

As it folloius from (|4.9|l the Lie algebra generated by functions (|4.8p 
is isomorphic to direct sum u(2) © u(2) © u(2) of three copies of u(2). 
The map 

J(z + ,z,v\v) :=0<f^r f*^) (4.15) 

is a Poisson map (momentum map) of (C 6 , { • , • }) into (u(2)*©u(2)*© 
u(2)*,{-,-} L p), i.e. 

{fo J, go J} = {f,g} LP o J, (4.16) 

uuhere noiu by { • , -}lp me denote sum of three copies of Lie Poisson 
bracket (|4~T3)| . 

Writting the Hamiltonians (|4.2|l - (|4.7|l in terms of the variables (|4.8|l 
me obtain on u(2)* © u(2)* © u(2)* the folloming three Hamiltonians 
in involution 

2 

R := a a 1 = r 1 + r 2 + r 3 (4-17) 

Ct 2 — 0-1 



# = (vi + V2 + + m + m) + J(n + r 2 + r 3 ) 2 + (4.18) 

+ -(ai - a 2 )(d 1 r 1 + d 2 r 2 + d 3 r 3 ) + ^(si + s 2 + s 3 ) 2 + 
+ -(ai + a 2 )(rfisi + d 2 s 2 + C/3S3) 
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G := F - (at + a 2 )H + ^(a 2 - ax)(sx + s 2 + s 3 )i2 + aia 2 £i- (4.19) 

- -(d + a 2 )(^Si + fi£js 2 + cfgS 3 )- 

- - (si + s 2 + s 3 )(disi + d 2 s 2 + d 3 s 3 ) = 

= (Vi + V2 + V3){difji + d 2 f)2 + d 3 fj 3 )+ 
+ {fji + m + fk)(di7i! + d 2 7] 2 + d 3 r} 3 )+ 

+ -(ax - a 2 )(d\rx + d\r 2 + dlr 3 ) + 
+ ^(n + r 2 + r 3 )(dxrx + d 2 r 2 + d 3 r 3 ) 



(4.20) 



and the Casimirs 



So = sx + s 2 + s 3 (4.21) 

5x = dxSx + d 2 s 2 + d 3 s 3 (4.22) 

5 2 = dj Sl + djs 2 + d 2 3 s 3 (4.23) 

uuhich are expressed by the Casimirs s±, s 2 , and S3. 
Since one has 

l -c k = sl = rl + (2\T lk \) 2 (4.24) 

for k = 1,2,3 we conclude from (|4.2ip - (|4.23p that if s k 7^ then the 
Hamiltonian hierarchy (|4.2 |l -(|4.7 |l is reduced to the three Hamiltonians 
(|4.17p - (|4.19p defined on the six dimensional symplectic leaves 

S si ,s 2 ,s3 := ^si x &s 2 x S 2 3 (4-25) 

mhich are the products of tiro dimensional spheres uuith radii Sx, s 2 , 
and S3 respectively. 

In order to simplify notation me mill use the polar coordinates r k , <p k 

on S sljS2)S3 defined by 



Vk = ^sl-rl^K (4.26) 

Note here that r k and <p k , —s k < < s k and — | < <ft k < f , mhere 
k = 1, 2, 3 and s k > 0, form canonical system of coordinates (Darboux 
coordinates) on E Sl)S2|S3 , i.e. 

{r k Ai} = hi- (4.27) 
From (|4.27p it folloms that the symplectic form on S SliS2iS3 is given by 
co Sl , S2 ,s 3 = drx A d(j)x + dr 2 A d(j) 2 + dr 3 A d(p 3 . (4.28) 
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In the polar coordinates the Hamiltonians l|4.17jl - l|4.19)l take the follom 
ing form 



H 

2 
1 



i^ (s 2_ r 2 )(s 2_ r 2 )cos(0i _ h)+ (4 29) 



+ \>J(4 ~ 4)(s 2 3 - r$) cos(fa - fa) - -{r\ + r\ + r 3 2 ) + 
+ \{s\ + si + s\) + l( n + r 2 + r 3 ) 2 + 
+ -(ai - a 2 )(dir a + d 2 r 2 + d 3 r 3 ) + ^(s 1 + s 2 + s 3 ) 2 + 
+ -(d + a 2 )(d x si + d 2 s 2 + d 3 s 3 ) 

C7 = i(d a + d 2 )^( s ?_ r 2)( s 2_ r 2 )cos(0i - 2 )+ (4.30) 
+ \{di + d 3 )^(sl-rl)(s 2 3 -r 2 3 ) cos(fa - 3 ) + 
+ ^2 + d 3 )^(s 2 - r 2 )(s 2 - r 2 ) cos(0 2 - 3 )- 
- ^(^i r i + ^2 + d 3 rf) + + d 2 s\ + c/ 3 s 2 ) + 

+ i(ai - a 2 ){d\r x + d 2 r 2 + d 2 r 3 ) + 
+ -(ri + r 2 + r 3 )(diri + <i 2 r 2 + <i 3 r 3 ) 

J R = r 1 +r 2 + r 3 . (4.31) 

The function i? generates on the symplectic leave S Sl)S2>S3 the follow 
ing Hamiltonian flow 

crf(r u r 2 , r 3 , fa, fa, fa) = (n, r 2 , r 3 , fa + <p,fa + <p, fa + fa. (4.32) 

Reducing the Hamiltonian system (Y} SljS2!S3 ,U) SuS2iS3 , H,G, R) to the 
level i? _1 (r) of the function R defined in (|4.31|l me obtain the reduced 
phase space R^ 1 {r) / {a^} uuith symplectic form oj given by 

w = dri A dfa + dr 2 A dfa (4.33) 

inhere fa := fa - fa and fa := fa - fa. 
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By introducing the notation 



/(ri,r 2 ):=-V( a ?-r?)( a 2-r§) (4.34) 
1 



fi(ri,r 2 ) := -yj (sj - r?)(s| - (r - n - r 2 ) 2 ) (4.35) 

/ 2 (ri,r 2 ) := 2 - r i)( s i ~ ( r ~ r i ~ r s) 2 ) ( 4 - 36 ) 

/i(n, r 2 ) := -j( r i +rl + (r-r 1 - r 2 ) 2 )+ (4.37) 

+ ^(ai - a 2 )(diri + d 2 r 2 + d 3 (r - n - r 2 ))+ 

+ ~(si + s 2 + s 3 ) 2 + ~(cn + a 2 )(disi + d 2 s 2 + tZ 3 s 3 ) 

(/(n, r 2 ) := -^(c/ir 2 + d 2 r 2 + d 3 (r - n - r 2 ) 2 )+ (4.38) 

+ ^(rfis 2 + d 2 s 2 + d 3 sl) + 

+ -(ai - a 2 )(d 2 ri + <i 2 r 2 + ^(t" - n - r 2 ))+ 
r 



u;e can uurite the Hamiltonians (|4.29|l and (|4.30|l as follouus 



# = f(r 1 ,r 2 ) cos(^i - ^2) + h(ri, r 2 ) cos(^i)+ 

+ Mri, r 2 ) cos(^ 2 ) + r 2 ) (4.39) 

G = (di + d 2 )f(r 1 , r 2 ) cos(^i - V> 2 ) + (di + r 2 ) cos(^i)+ 

+ (c?2 + d 3 )f 2 (ri,r 2 ) cos(^ 2 ) + £/(n, r 2 ). (4.40) 
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From the above formulas rue find that the Hamilton equations for 
the Hamiltonian (|4.39p in the canonical coordinates 7*1, r 2 , 0i, ip 2 as- 
sume the follouuing form 

^ = {^1, H} = -|£ cos(Vi - V2) - ^ cos(Vi)- (4.41) 
ari ari 

cos(t/>2) - 

r\ or 1 

2 = {0 2 , H} = -^- cos(Vi - V2) - ^ cos(Vi)- (4.42) 
ar 2 ar 2 

_ ^ cos(0 2 ) - — 
<9r 2 <9r 2 

ri = {ri,i?} = -/(ri,r 2 )sin('0i - ^2) - fi(n, r 2 ) sin^i) (4.43) 

f 2 = {r 2 , if} = /(ri, r 2 ) sin (-01 - 02 ) - Mri, r 2 ) sin(0 2 ) (4.44) 

Equations (|4.39|| - (|4.40|| define 0i and 2 as an implicit function of 
ri and r 2 . In order to solve equations (|4.4ip - (|4.44p me apply the gen- 
erating function method. Thus me introduce nem variables 7 and r 
canonically conjugated to the integrals of motion G and H, i.e. 

u = -dfyidn + 4> 2 dr 2 ) = d(>ydG + rdH). (4.45) 

From (|4.45p it folloms that there exists a locally defined function $(r 1; r 2 , G, 
such that 

d$ = ^) x dr x + il) 2 dr 2 + ^dG + rdH. (4.46) 
The equality (|4.46p gives the relationship betmeen nem canonical co 
ordinates (7, r, G, H) and the old canonical coordinates (r 1; r 2 , 1; r/> 2 ) 
given that me obtain the generating function $. To this end let us 
note that consistency condition 

^1 = ^ (4.47) 

ar 2 ari 

for the equations 

<9r 2 

folloms from {iT, C7} = 0. The proof of this fact is a consequence of 
the formula on the derivative of the implicit function 



01 (4.48) 

02 (4.49) 



<902 <902 
\ dr\ dr2 J 



dG dG\ - 1 /dG dG 



-(M it [k M\- <4 - 5o) 

,dipi (902/ \dri dr2, 
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Thus me can define $ for (r 1 ,r 2 ) G [— Si,Si] x [— s 2 ,s 2 ] by 

i 

$(ri,r 2 ,G,H) = - J (if>i(sri, sr 2 ) n + i/j 2 (sr 1 , sr 2 ) r 2 )ds. (4.51) 
o 

Recall that the dependence of ipi and ip 2 on r 1; r 2 , G, H is given in 
the implicit may by l|4.39jH|4.40j> . 

Nom submitting $ defined by (|4.51|l into 

m = T < 4S3 > 

me find the nem coordinates 7 and r. 

Since the variables (7, r, G, if) are angle-action coordinates for the 
Hamiltonian system defined by H me obtain 

r(t) = t + t (4.54) 

7(t) = const (4.55) 

G(t) = const (4.56) 

H(t) = const (4.57) 

The time dependence of the variables ri(t) and r 2 (t) can be found 
in the implicit may by l|4.52) - <|4.53) . Next me obtain the time depen 
dence of ipi(t) and ijj 2 (t) from (|4.48p (|4.49p . 

Summing up, me have solved Hamiltonian system given by Hamil 
tonian (|4.2p describing nonlinear interaction of six maves. The solution 
mas obtained in quadratures, but due to technical difficulties, explicit 
form of solution mould be too complicated to present. 

5. Solution in (2 + 2) dimensional case 

In this section me consider a special situation of (2 + 3)-dimensionaI 
case, mhen S3 = 0. From (|4.24p me see that S3 = implies that r 3 = 
and I ^73 1 = 0. Since S3 is Casimir, this anzatz is consistent mith the 
evolution mith respect to all Hamiltonians (|3.3p . 

Therefore me mill solve equation (|3.23|l in the (2 + 2)-dimensionaI 
case taking instead of 6 dimensional symplectic leaves (|4.25p , the 4- 
dimensional symplectic leaves 

S SljS2 := x S^ 2 . (5.1) 

mith the symplectic form given by 

(jJ Si ,s 2 = dri A d(j)\ + dr 2 A d(f> 2 . (5.2) 
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The Hamiltonians l|4.17jl - l|4.19)l assume nouu the form 

H = \^{s\-rl){sl-rl)cos{<j> x - 2 ) - \{r\ + r 2 2 )+ (5.3) 
+ \(4 + si) + \(n + r 2 ) 2 + - a 2 )(d iri + d 2 r 2 ) + 
+ ^(si + s 2 ) 2 + i(ai + a 2 ){d 1 si + d 2 s 2 ) 



G=\{d 1 + d2)^/(s?-r?)( S l-r|) co8(^i - (5.4) 

- ^(^l^i + c^r 2 ) + ^(^lSi + 4^2) + -(at - a 2 )(^i + d 2 2 r 2 ) + 

+ 7^1 + r 2 )(diri + d 2 r 2 ) 
R = n + r 2 . (5.5) 

Note that in this case no longer rue need G as an additional in- 
tegral of motion, since in (2 + 2)-dimensionaI case the integrals of 
motion H and R are sufficient to integrate the system. Similarly 
to the (2 + 3) dimensional case, me reduce the Hamiltonian system 
(E SljS2 , U) SliS2 , H, R) to the level set _R _1 (r) and obtain the Hamiltonian 
on the reduced phase space i? _1 (r)/{cr^}^ eR 

H = ^ v /( S ?-r2)( s 2_( r _ ri )2 )cos(V , l) + W2(ri)) (5 . 6) 

inhere 

w 2 (n) := --{rl + (r - n) 2 ) + ^(s\ + s 2 2 )+ (5.7) 
+ ^r 2 + ^(ai - 02) (din + d 2 (r - r x ))+ 
+ ~0i + s 2 ) 2 + ~(ai + a 2 )(disi + d 2 s 2 ). 

The symplectic form uj on i2 -1 (r)/{a - ?}^ £ R is given by 

w = dri A dV>i, (5.8) 

mhere ipi := 0i — 2 - The flouj {a^j^R defined by the integral of 
motion R has the form 

of (ri, r 2 , 0i, 2 ) = (ri, r 2 , fa + 0, 2 + 4>)- (5.9) 
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The equations of motion (|4.41|l - (|4.44|l in this case reduces the fol 
louring ones 

i) X = {if> u H} = n-~r+ i(oi + o 2 )(di - d 3 )- (5.10) 

-r 1 (s^-(r-r 1 ) 2 ) + (s?-rf)(r-r 1 ) 

— cos tpi 



fi = {ri,#} = - r ?)( s 2 - (r-r^^sin^!. (5.11) 

In order to solve the equation (|5.1 1)1 me use (|5.6p and Pythagorean 
identity to obtain the relation 

A(H - w(n)) 2 + 4(r x ) 2 = {si - r 2 )( S 2 - (r - n) 2 ). (5.12) 

The solution of (|5.12p is in the form of elliptic integral of the first kind 

* = ± (5.13) 

mhere 

tw 4 (ri) := (sj - r 2 )(s 2 - (r - n) 2 ) - A{H - w 2 {r x )f (5.14) 

is a polynomial of fourth order. Thus r\ is an elliptic function of the 
parameter t. Subsequently, from (|5.6p me obtain 

ipi = arccos — = 2 ^ 1 ^ (5.15) 

V (4 - r l)(4 - (r - n) 2 ) 

In the analogous may me can find solution of the Hamiltonian sys 
tern given by G or, more generally, for any Hamiltonian or H 5 \. 

6. Physical interpretation of the (2+3) mode Hamiltonian 

In order to elucidate the optical interpretation of the Hamiltonian 
(|4.2|l let us express it in the coordinates z%, z 2 , z% and v%, v 2 , v%. 



2 

H = a\dx \z\\ + a\d 2 


W 


2 + airf 3 


1^3 




+ a 2 d 2 \v 2 2 + a 2 d 3 


1^3 


2 + 






l/i 1 A 1 ,4 

+ -(ki| + N + 




| 4 +h| 4 


+ 


1^2 1 + |^3] ) + 


1 1 - 1 1 2 . 1 |2 


N 


2,1 |2 




2,1 |2 1 |2 , 1 i2 1 |2 , 1 i2 1 i2 , 
+ Fl| P2I + Pl| |f 3 | + M |f 3| + 


1 1 |2 1 |2 , 1 |2 
+ |^l | tll| + \z 2 \ 


1^2 


2,1 i2 


|«3 


| 2 + 



+ ViV 2 Z 2 Zi + V 2 ViZ±Z 2 + V^Z^Zl + V 3 V 1^1 ^3 + V 2 V 3 Z 3 Z 2 + V 3 V 2 Z 2 Z 3 , 
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Assuming that these coordinates describe separate modes of the six- 
luave interacting nonlinearly through the nonlinear dielectric medium, 
rue get the following optical interpretation of the particular terms: 

i) the quadratic terms \zi\ 2 and \vi\ 2 constitute free Hamiltonian Ho, 
i.e. they describe free energy of the light, inhere a^dj is propor 
tional to the corresponding modes frequency; 

ii) the terms |zj| 4 and \vi\ 4 are responsible for the Kerr effect, i.e. 
third order of nonlinearity of polarisation of a medium causing 
intensity dependent phase shift, see e.g. [WM95J; 

iii) the terms \zi\ 2 \zj\ 2 , \zi\ 2 \vj\ 2 , \vi\ 2 \vj\ 2 introduce Kerr-Iike effect, 
i.e. phase shift of i th mode depending on intensity of j th mode; 

iv) the other terms describe the conversion between the modes, 
e.g. the term v\ViZiZ\ describes the process of absorption by 
the medium of certain amount of light in modes Z\ and t>2 with 
simultaneous emission of light in modes V\ and z-i- 
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